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Abstract 

The evolving galaxy is considered as a system of baryonic fragments em- 
bedded into the static dark nonbaryonic (DH) and baryonic (BH) halo and 
subjected to gravitational and viscous interactions. Although the chemical 
evolution of each separate fragment is treated in the frame of one - zone close 
box model with instantaneous recycling, its star formation (SF) activity is a 
function of mean local gas density and, therefore, is strongly influenced by 
other interacting fragments. In spite of its simplicity this model provides a 
realistic description of the process of galaxy formation and evolution over the 
Hubble timescale. 

1 Introduction 



Recent advances in extragalactic astrophysics show the close link of a disk galaxy 
dynamical evolution and its chemical and photometric behaviour over the Hubble 
timescale. In spite of remarkable succes of the modern theory of galaxy chemical 
evolution in explaining the properties of evolving galaxies (M) its serious shortcom- 
ings concern the multiparameter character and practical neglecting of dynamical 
effects. The inclusion of simplified dynamic into the chemical network (||) and vice 
versa the inclusion of simplified chemical scheme into the sophisticated 3D hydro- 
dynamical code (0, HJ) gives very promising results and allows to avoid a formal 
approach typical to standard theory. 

In this paper the interplay between a disk galaxy dynamical evolution and its 
chemical behaviour is studied in a frame of a simplified model which provides a 
realistic description of the process of galaxy formation and evolution over the cos- 
mological timescale. 
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2 Initial conditions 



The evolving galaxy is treated as a system of baryonic fragments embedded into 
the extended halo composed of dark nonbaryonic and baryonic matter. The halo 
is modelled as a static structure with dark (DH) and diluted baryonic (BH) halo 
components having Plummer - type density profiles (0): 



pBH{r) = T 



M BH b% H 



3 



7T 



and 

Pdmh{j) 



Mdmh ^dmh 



where 



M DMH = 10 12 M , b DMH = 25 kpc, 



M BH = 10 11 M Q , b BH = 15 kpc. 

The dense baryonic matter (future galaxy disk and bulge) of total mass of M gas = 
10 11 M Q is assumed to be distributed among N = 2109 particles - fragments. The 
single particle density profile is also assumed to be a Plummer one. Its mass is 
taken to be m 8 = M gas /N and radius hi = 1 kpc. Initially all particles are smoothly 
placed inside the sphere of radius R gas = 50 kpc and are involved into the Hubble 
flow (Hq = 75 km/s/Mpc) and solid - body rotation around z axis. The initial 
motion of this system is described as: 

V(ar, y, z) = [ft{x, y, z) x r] + H ■ r + DV(ar, y, z), 
where il(x,y,z) = (0,0,1) • Vt cir is an angular velocity of the rotating sphere, 

^cir ^cir / -Rgas and 



Vrir — 4 / G • 



M gas + Mdmh + M BH 



cir p 
It r , 



^gas 

The components DV X , DV y , DV Z of the random velocity DV vector are assumed 
to be initially randomly distributed within an interval -j- 10 km/s. 



3 N — body code 

The dynamical evolution of baryonic matter fragments which are subjected to grav- 
itational influences of DM-baryonic halo and interfragment interactions is followed 
by means of effective N - body integrator with individual time step. The dynamics 
of such N - body system is described by following equations: 

dUi/dt = Vj, 

(1) 

k dVi/dt = Ai(R,V). 
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The acceleration of i - th particle Aj is defined as a sum of three components. 



A, = A[ NT + Af XT + AY ISC , (2) 

where the first term A^ 7 accounts for gravitational interactions between frag- 
ments. The second one Af XT is defined as an external gravitational acceleration 
caused by the DM and baryonic halo. The last term A^ 180 corresponds to the 
viscous decceleration of fragment when passing through the baryonic halo. 

The gravitational interaction between fragments is defined as the interaction of 
N Plummer profile elements: 



v 



A{ NT = -G- V ^ g-.Ri,. (3) 

Here h{j = (hi + hj )/2 and R^- = Rj — R^. 

Accounting for that the halo DM and baryonic components are also Plummer 
spheres the second term becomes: 

Af XT = -G ■ ( Mdmh 3 + MBH 3 ) " ^ (4) 

The form of the last term A\ ISC will be discussed in the next subsection. 
The characteristic time step 5ti in the integration procedure for each particle is 
defined as: 



5tj = Const ■ min[ 
j 



R-zj I I R-ij 



(5) 



where the Vy = V, — Vj and the A^ = Aj — Aj. 



Here the Const is a numerical parameter equal to Const = 10~ 2 that provides 
a nice momentum and energy conservation over the integration interval of about 15 
Gyr. For example, in the conservative case (when viscosity of the system is set equal 
to 0), the final total error in the energy equation is less than 1%. 



4 Viscosity model 

The viscosity term AY ISC is artificially introduced into the model so as to match the 
results of more sophisticated SPH approach on dynamical evolution of disk galaxies. 
The best fitness of results of this simplified approach with SPH modelling data (see 
e.g. is achived when the momentum exchange between the baryonic halo and 
moving particles is modelled by the following expression: 

AJ ISC = —k-Vj- _LYi_L . PI"*! . . (6) 

1 E? n rrr>9 as I ™ Star \ > 

itylSC PVISC ™>i + 

Here pvisc and Rvisc are numerical parameters set equal to pvisc = 0.2 cm -3 
and Rvisc = 10 kpc. The vector Vj is a particle velocity vector. It is to be noted 
that single particle is assumed to have a total mass which doesn't change with time 
and is defined as sum mi = mf as + mf ar . But masses of its gas and star components 
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mf as and mf ar are variable values and are defined by the temporal evolutionary 
status of the given fragment. Initially mf ar = 0. The results of fitness show that for 
z component of viscosity term k — 1. In the galaxy plane where it is necessary to 
account for baryonic halo and baryonic fragments partial corrotation the dynamical 
friction is decreased and for x, y component of viscosity term k reduced to the value 
0.15. 



5 Density definition 

In the frame of the multifragmented model the definition of local gas density is 
introduced in the SPH manner, e.g. local gas density depends on the total mass of 
matter contained in the sphere of radius Hi around the i - th particle. For each 
i - th particle the value of its smoothing radius Hi is chosed (using the quicksort 
algorithm) requiring that the volume of such radius compraises Nb = 21 nearest 
particles (i.e. ~ 1 % of total number of particles N). Therefore, the total mass Mj 
and density of gas pi inside this sphere are defined as 

Mi = £ Amff, * = J^, (7) 



where Ara? os is defined as 



if | R y \>(Hi + hj) => Amff = 0, 
if | Kij \< (Hi - hj) => Amfj 8 = mf% 

else Amff = mf s ■ Eilhzl^iil , 

6 Star formation and SN explosions 

A forming disk galaxy is modelled as a system of interacting fragments (named as 
particles) embedded into the extended halo. Each particle is composed of gas and 
stellar components and its total mass is defined as m, = m\ as + mf tar . Initially all 
particles are purelly gaseous and, therefore, initially mf ar = 0. To follow a particle 
star formation (SF) activity a special timemark t b ^ 9SF is introduced which initially 
is set equal to t b f gSF = 0. The particles eligible to star formation events are chosen 
as particles which still have a sufficient amount of the gas component and their 
densities exceed some critical value p m inSF during some fixed time interval At^ (of 
order of free-fall time): 

Pi > PminSF, 

U - t begSF > AtsF, (8) 

k mf as /(m g i as + mf ar ) > 1CT 4 . 
Here At SF = 50 Myr, and 

PminSF — 0.01 cm 3 (this last value is not crucial and 

is only limiting one). 
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If the particle was subjected to SF activity the parameter t b ^ 9 F is set equal to 



^begSF = f.^ aric [ m ^ tar anc l mf as are redefined as 

r^star r /i T)\ ^9 as i ^star 

m i .— e ■ [L — -ti) ■ vn i + m i , 

mf as := (1 — e ■ (1 — R)) ■ m\ as . 
Here e is a SF efficiency which is defined as 



(9) 



e = a.-^.(l-exp(-— ))• (10) 

PSF Pi 

Therefore, 

■ t Pi n ^ Pi 
ij > => e — > a ■ , 

PSF PSF 

if ^£^ =>e^a. 

Pi 

To match available observational data on star formation efficiency (see e.g. ||) 
parameters a and p$F are set equal to a = 0.5 and psf = 10cm -3 . The chemical 
evolution of each separate fragment is treated in the frame of one - zone close box 
model with instantaneous recycling. In the frame of this approach for the returned 
fraction of gas from evolved stars is taken a standard value R = 0.25. 

Following the instantenuous recycling approximation it is assumed that after 
each SF and SN explosions the heavy element enriched gas is returned to the system 
and mixed with old (heavy element deficient) gas. After each act of SF and SN 
explosions the value of heavy element abundances of gas in particle is upgraded 
according to: 

e-R-AZ 

Zi := Zi + ; r. 11 

1 1 - e • (1 - R) K J 

The value AZ = 0.01 is used as an average value for all values of Z = 0.001^0.04 
(see 0). Initially Zi = 0.0 in all particles. 



7 Conclusion 



The proposed simple model provides the self - consistent picture of the process of 
galaxy formation, its dynamical and chemical evolution is in agreement with the 
results of more sophisticated approaches (see e.g. [0, H 1, 0]). 



• The rapidly rotating protogalaxy finally formed a three - component system 
resembling a typical spiral galaxies: a thin disk and spheroidal component 
made of gas and/or stars and a dark matter halo. 

Fig.l and Fig. 2 show respectively the star formation rate and the galaxy total 
stellar mass as a function of time. Fig. 3 shows the cylindrical distribution of 
stellar (upper curve) and gaseous (lower one) components of the final model 
disk galaxy as a function of a distance from the galaxy center in the galactic 
plane. 
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The total star formation rate (SFR) is a succesion of short bursts which doesn't 
exceed 28 M Q /yr. During first 2 Gyr of evolution only about 20% of total 
galaxy mass is transforms into the stars. The SFR gradually decreases, during 
the further evolution, to the value of about 2 M /yr typical for our own 
Galaxy. The final total stellar and gas mass of the model galaxy disk are 
about 92% and 8%. All these data as well as surface densities distributions of 
stellar and gaseous components (see Fig.4) are in nice agreement with present 
date observational data (|| U). 

• The disk component posesses a typical spiral galaxy rotation curve and the 
distribution of radial and V z - th velocities of baryonic particles clearly show 
the presence of the central bulge (see Fig. 6 and Fig. 7). 

• The metallicities and the global metallicity gradient resemble distributions 
observed in our own Galaxy (Fig. 5). The averaged observed value of global 
metallicity Z/Z@(r) (see 0]) is shown in this Fig. as a solid line. 

Acknowledgements: Peter Berczik would like to acknowledge the American 
Astronomical Society for financial support of this work under International Small 
Research Grant. 

8 Figures 

Figure 1: The variation of total star formation rate of forming disk galaxy with 
time. 

Figure 2: The growth of galaxy stellar mass with time. 

Figure 3: The cylindrycal distribution of masses of stellar (upper curve) and gaseous 
(lower one) components as a function of distance from galaxy center in a galactic 
plane. 

Figure 4: The surface density radial distributions of stellar and gaseous components 
(stellar component is shown by filled stars, gaseous one by astericks). Theoretical 
distributions (not scaled) of surface density for radial exponential scale lengthes 2.0 
and 3.0 kpc are shown below by lines. 

Figure 5: The radial distribution of heavy element abundances (averaged observed 
distribution of z for our Galaxy is shown by solid line). 

Figure 6: The final distribution of V z - th velocities of baryonic gas-stellar particles. 
Figure 7: The final galaxy rotation curve. 
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